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Abstract. We present a detailed study of the scaling behavior of correlations functions and AC susceptibil- 
ity relaxations in the aging regime in three dimensional spin glasses. The agreement between simulations 
and experiments is excellent confirming the validity of the full aging scenario with logarithmic corrections 
which manifest as weak sub-aging effects. 
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1 Introduction 

There is a great interest in the understanding of dynam- 
ical effects in spin glasses. These include magnetization 
relaxation, aging and temperature change protocols 

The study of these effects may clarify the nature of 
spatial effects and coarsening phenomena in spin glasses, 
an issue which remains still poorly understood. 

In this paper we present a detailed study of magne- 
tization relaxation phenomena and aging effects in three 
dimensional Edwards- Anderson spin glasses. Our primary 
goal is to check the validity of the full t/tw scaling behav- 
ior in correlation functions as well as identifying possible 
sources of corrections to that behavior by comparing to 
experimental data. Dynamical experiments in spin glasses 
include magnetic relaxation and AC measurements. Here 
we will focus our attention on correlation function and 
AC susceptibility relaxations. The advantage of studying 
correlations is that these are easy to evaluate numerically 
being also tightly related to thermoremanent magnetiza- 
tion relaxation experiments. On the other hand, AC relax- 
ations can be directly compared to experimental results 
and, to the best of our knowledge, no results appeared on 
this point in the literature. 

There exist many works on simulations in the litera- 
ture studying correlations or remanent magneti- 
zation relaxations and this part of the topic that we in- 
vestigate here is certainly not new. What has never been 
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considered in detail in the past and merits further investi- 
gation is the explicit comparison between simulations and 
experiments. Having the experimental results in mind we 
have tried to apply the same scaling plots used by the 
experimentalists to our numerical data. This may serve 
as a valuable guide to better understand what properties 
are generic to spin glasses and what coarsening scenario 
accounts for the collected experimental data. 

Magnetic relaxation (or correlation function) and AC 
experiments give equivalent information, the advantage of 
using AC experiments is that they constitute a very sensi- 
tive tool to detect dissipative processes. When measuring 
correlations the external timescale tw is fixed by the time 
elapsed after quenching below Tc while in AC experiments 
the external timescale is fixed by the inverse of the fre- 
quency of the AC field. In a full scaling scenario, in the 
first class of experiments the relevant scaling variable is 
t/tw while in the second class it is ut. In what follows we 
check the validity of this simple scaling behavior identify- 
ing possible sources of corrections. 



2 The model and the observables 

The Edwards- Anderson model ||l^ is defined by the fol- 
lowing Hamiltonian 



n 



i=l 



(1) 



where the indices i, j run from 1 to V, the ai are Ising spins 
{(7i = ±1) and the pairs identify nearest neighbors 
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in a three dimensional lattice. The exchange couphngs Jy 
are taken from a random distribution. The simplest choice 
is a Gaussian distribution with zero average and finite 
variance, 



(2) 



This model displays a spin glass transition at finite 
temperature Tc ~ 0.95 [|T^,|T2|. A Monte Carlo step cor- 
responds to a sweep over V randomly chosen spins of the 
lattice. Monte Carlo simulations of (|l|) use random updat- 
ing of the spins with the Metropolis algorithm. Dynamical 
experiments use very large lattices (typical sizes are in the 
range L = 20 — 100) with negligible finite-size effects for 
the largest sizes {L = 64 for magnetization relaxation ex- 
periments and L = 100 for AC experiments). Correlation 
function simulations have been done on a special purpose 
machine APEIOO ||l^ for sizes 64'^ and averaged over 10 
samples. AC experiments were done for a single sample on 
a Linux cluster of PC's for size L = 100. 

Relaxation measurements are done applying a uniform 
magnetic field and measuring the decay of the thermore- 
manent magnetization (hereafter denoted by TRM), or 
equivalently, the growth of the zero-field cooled (ZFC) 
magnetization. The typical experiment consists in the fol- 
lowing. A sample is fastly quenched below the spin glass 
transition temperature for a time tw (i.e. the waiting time). 
Then a uniform small magnetic field h is applied and the 
growth of the magnetization measured, 



XZFc{tw,tw + t) 



1 ^ . 



Vh 



t) 



(3) 



Another quantity of interest related to the magnetiza- 
tion which can be numerically investigated is the two-time 
correlation defined by 



C{ty^ , -\- t) 



i—1 



(4) 



The interest of studying correlations instead of zero- 
field cooled magnetizations is that they yield the same 
dynamical information. Indeed in the stationary regime 
they are related through the fluctuation-dissipation theo- 
rem (FDT) 

XzFc{t) = ■ (5) 

In AC experiments an oscillating magnetic field h{t) = 
ho cos(27rci;t) of frequency uj = ^, where P is the period, 
is applied to the system and the magnetization measured 
as a function of time 



M (t) = Mo cos(27rwt - 



(6) 



where Mq is the intensity of the magnetization and is the 
dephasing between the magnetization and the field. The 
origin of the dephasing is dissipation in the system which 
prevents the magnetization to follow the oscillations of the 



magnetic field. From the magnetization one can obtain the 
in-phase and out-of-phase susceptibilities defined as 

, _ Mq cos(0) _ 2 /(f M{t) cos{2TTUjt)dt 

ho ho 
„ Afosin(0) 2J^'' M{t)smi2TTu;t)dt 

^ ho ■ 

The dephasing (p measures the rate of dissipation in the 
system and is given by 



tan(0) - ^ 
X 



(9) 



In mimerical simulations the in-phase and out-of-phase 
susceptibilities are computed by averaging the right-hand 
side in Eqs.(0),(^) over several periods P — ^ - This means 
a very large measurement time for low frequencies for both 
experiments and simulations. 

In the numerical simulations (both in DC or AC exper- 
iments) the intensity of the probing fields cannot be arbi- 
trarily small because of the weakness of the signal in com- 
parison to other source of fluctuations such as finite-size 
effects (which induce finite- volume statistical fluctuations 
for extensive quantities, like the susceptibility). Conse- 
quently, the intensities of the probing magnetic fields used 
in numerical simulations are much larger than the corre- 
sponding experimental ones (between 50 and 500 times 
larger). As we will comment later, we do not believe that 
this leads to conflicting results between simulations and 
experiments. As soon as one checks that measurements 
are done within the linear response regime then the in- 
tensity of the probing field should not be crucial. Actu- 
ally the values of the intensity of the fields usually em- 
ployed in numerical simulations of TRM experiments are 
well known to satisfy linear response ||lj,|l5|. Note that, 
in general, similar difficulties are encountered when ana- 
lyzing data in both numerical simulations and real exper- 
iments, the main difference is the absolute magnitude of 
the time scales one can explore in the two cases (up to mi- 
croseconds in simulations and between seconds and days 
in experiments). 



3 TRM and correlation function relaxations 

In order to study time scaling in a wide times range (spe- 
cially in the aging t ^ tw regime) experimentalists have 
measured the decay of the TRM which is strictly re- 
lated to the zero-field-coolcd one through Mzfc = Afpc ~ 
A/trm, where the field-cooled magnetization is practically 
constant in the glassy phase. 

From the numerical simulations point of view the best 
quantity one can look at for checking time scaling is the 
autocorrelation function C {tw ,tw -\- t). It has much less 
fiuctuations than any response to an external field. More- 
over in the quasi-equilibrium regime where the fluctua- 
tion-dissipation theorem (FDT) holds, it gives exact in- 
formation on the zero-field-cooled susceptibility via 



+ t) 



1 C(iui J tiu + t) 
f 



(10) 
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Fig. 1. Aging at T = 0.5 with and without an external mag- 
netic field after time for a 64^ system. Waiting times (t™) 
range from 10^ (leftmost curves) to 10^ (rightmost curves). 



In the aging regime the connection between correlation 
and susceptibihty is less trivial. However in the limit of 
small perturbing field and large times, where a general- 
ization of the FDT seems to hold |l^,|r^,|l5|, a relation 
between correlation and susceptibility can still be estab- 
lished. Then, in general, concerning time scaling, we can 
safely assume that the autocorrelation functions decay as 
the TRM do. 

In this section wc try to understand which is the best 
scaling for the C (i„ , -I- data. The measurements have 
been taken on 10 samples of a 64'^ system, at a tempera- 
ture T = 0.5, with waiting times ranging up to — 10^ 
and measuring times up to t = 10®. We have considered 
two different experimental situations, that is with or with- 
out an external magnetic field after time tw (the evolution 
up to time being always with no field) , in order to check 
whether such a small perturbation may change the dy- 
namical scaling. The external field intensity /i = 0.1 has 
been chosen such that the system is in the linear response 
regime. If the magnetic field would be applied during all 
the experiment we do not expect any sensible difference 
with the ft, = case. 

In Fig. we show the correlation functions data, with 
and without the external magnetic field. As the time goes 
on the effect of the magnetic field seems to accumulate 
and the differences become larger. Note however that, for 
any given waiting time t^,, the correlation curve presents 
the two well known regimes [^|l^: the quasi-equilibrium 
one (t < t^) and the aging one (F > t^)- 

Because we are mainly interested in the scaling in the 
aging regime, we have tried, as the simplest analysis, to 
collapse the t > data using t/t^ as the scaling variable. 
The results are shown in Fig. g and they clearly show 
that a value for /i smaller than 1 is needed in order to 
collapse the data in the large times limit. Note also that 
the presence of an external field seems to decrease sensibly 
the value of /i. The errors on the estimation of /i are of 
the order of 10~^. 



This numerical result may suggest the presence of a 
sub-aging regime in the EA model |l9| and it could be 
interpreted as one more similarity with real spin glasses, 
where /i = 0.97 j^]. However a more careful analysis shows 
that the correlation functions are perfectly compatible 
with a full aging, that is t/tw, scaling. In Fig. ^ we present 
the results of such an analysis, which has been done fol- 
lowing the one performed on experimental TRM data in 
reference 

Few comments are in order. The best values for the A 
and a parameters seem to be similar to the experimen- 
tal ones {A ~ 0.1 and a = 0.02). However, because of 
the lack of a quantitative criterion for data collapsing, the 
best collapse is very often subjective. In this case we have 
found that, in order to obtain a good data collapse, the 
parameter must be fixed to 1 or very close to it. On the 
other hand, the A and a parameters are strongly corre- 
lated (with a correlation coefficient close to -1) and they 
can be changed by a quite large amount without affecting 
the data collapse. Then the errors on these parameters are 
large. In Fig. y we show the collapse for parameters values 
being more or less in the center of the confidence region. 
We have also tried to collapse both sets of data {h ~ and 
h — 0.1) with the same parameters, but it was impossible 
to obtain any reasonable data collapse. 

In Fig. ^ we use the scaling variable log(i -I- i^,) — 
log{t^), even if t/tw would be the most natural one when 
full aging holds. Our choice is dictated by the need for a 
comparison with the collapse of experimental data shown 
in Fig.S.b of reference ||]. There the scaling variable [{t + 
twY^'^ — ii,~^]/(l — /i) is used, which tends to \og{t + 
tw) — log(tm) in the fi 1 limit. It is well known that 
the goodness of a data collapse may depend on the scales 
chosen for presenting the data. In the present case, in the 
scaling variable \og{t + 1^^) — log{tm) we have better col- 
lapses that in the variable t/tyj, because large times are 
"compressed". Note however that both scaling variables 
give very good collapses of our data, the same being true 
for the experimental data . Anyhow it is worth to note 
that the use of the scaling variable t/t^ for checking full 
aging and [{t + t^Y^^^ — ^^7^1/(1 — ^) for checking sub- 
aging makes the life harder to the full aging scenario. This 
without considering the fact that there could be additional 
logarithmic corrections to the full t/tw scaling p9[ |. 

We conclude that it is not easy to obtain precise quan- 
titative information in order to distinguish the full aging 
from the sub-aging scenario. Moreover the presence of an 
external magnetic field, which on a first simple analysis 
seems to change the scaling, is in fact irrelevant for the 
scaling and it only changes a little bit the fitting parame- 
ters. 



3.1 A small note on the ZFC susceptibility scaling 

The scaling of the zero-field cooled susceptibility has been 
already studied in the past, both directly or via the 
fiuctuation-dissipation relation which links it to the cor- 
relation functions scaling |l^. Here we do not repeat this 
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Fig. 2. Best scaling in the aging regime {t > t^) for the data presented in Fig. |l| 
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Fig. 3. Full aging {/j, — 1) data collapse for the correlations functions presented in Fig. Note that the use of a logarithmic 
scale improves the quality of data collapsing. 



kind of analysis. We simply would like to present some 
new data regarding the ZFC susceptibility scaling. 

Indeed very recently Bernard! et al. have proposed the 
following scaling for the ZFC susceptibility |2l[] 

+t) = x{R{t^),m) , (11) 

where both length scales grow in an algebraic fashion 



R{t) oc L{t) cx t"^ 



(12) 



with an exponent linear in the temperature, like for the 
off-equilibrium correlation length |^,^,||]. We found that 
numerical data obtained with very large simulations are 
not compatible with the proposed scaling (see Fig. ^). 

Our numerical experiments are performed in the fol- 
lowing way. For any given temperature Ti (which takes 
3 values in our case Ti — 0.3,0.5,0.7), we start the sim- 
ulation from a random configuration and we let evolve 
the system at temperature Ti for a number of MCS t^, 
such that t = A where A is a constant that we fbced to 
A = 10^-^. At this time we switch on a small perturbing 
field, we move the temperature to T2 = 0.5 and then we 
measure the response of the system (ZFC susceptibility). 
From this kind of experiment we have obtained many in- 
teresting information on the spin glass low ternperature 
dynamics which have been published elsewhere |24|. Here 




Fig. 4. ZFC susceptibilities after a particular experiment (see 
text). If the scaling in Eq.(|ll[) would be correct, data in the 
figure should collapse. 



we used again this kind of experiment in order to verify 



the scaling ( |ll[ ) proposed in Ref. [21 



By construction we have that in all the three experi- 
ments the length R{t^) is the same and the L{t) is related 
to t by the same law, because after time we always make 
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evolve the system at the same temperature. Then if the 
scaling (|ll|) would hold, we should find a good data col- 
lapse in Fig. |[ which is not true. This behavior can be 
explained by observing that after the thermalization time 
tw{Ti) the three experimental situations are not identical: 
they have developed a similar correlation length, neverthe- 
less the actual configuration is different and the response 
to an external perturbation differs. 

In Ref. [21 the authors find a perfect agreement to 
the scaling (11). However it must be noted that they use 
temperatures in a small range, T e [0.5,0.7], which cor- 
responds to the two uppermost curves in Fig. B which 
indeed almost coincide. Violation to the scaling (pi] ) can 
be seen only at lower temperatures, which where not used 
in Ref. I2I. 



4 AC susceptibility relaxations 

Let us briefly remind how these experiments are usually 
done. The system is quenched to a low temperature (rang- 
ing between 0.6 and 0.9 times the value of Tg) and the AC 
susceptibility is recorded. Typical values of the frequency 
of the AC field are between 0.1 and 1 Hertz. The amplitude 
of the probing AC field is of order 10~^ to 10~^ Oersteds 
deep inside the linear response regime. The AC suscepti- 
bility is then recorded as a function of time and relaxation 
is observed on time scales of order ut ^ 1000 — 5000 corre- 
sponding to several thousands of periods of the AC field. 
Measurements are then obtained averaging over several 
cycles of the AC field in order to obtain x' and x" with 
enough accuracy (10 cycles is a typical value). 

Having in mind the experimental setup we have done 
the following experiment. Starting from a random config- 
uration we have measured the AC susceptibility as a func- 
tion of time for different frequencies of the AC magnetic 
field. The frequencies of the field are defined a,s u; — -p 
where P is the period of the oscillating field in Monte 
Carlo steps. The results for the in-nhase and out-of-phase 
susceptibilities are shown in Fig. |5| at temperature T — 
0.6 ~ 0.63 Tc and field periods P = 50,100,500,1000. 
Each point in Eqs.(^ and (||) is obtained by averaging over 
10 periods of the field. The curves can be very well fitted to 
a power law decays of the type xi^^ = xi^i 00) -I- At~'^ 
but the exponents depend very much on the frequency. For 
the largest frequencies the exponent in both susceptibili- 
ties are compatible with a value smaller than 0.1 although 
it is difficult to establish its precise value. A similar diffi- 
culty is found in laboratory experiments. 

To make evident the ujt scaling for the AC experiment 
we show in Fig. ^ the results of Fig. ^ plotted vs. ujt. The 
values of the susceptibilities in the vertical scale have been 
shifted by an arbitrary quantity to make them coincide. 
This procedure is exactly the same as done in experiments 
in the figure 2 of reference Hshowing the same qualita- 
tively results. Note from FigTH that the relaxation of the 
AC susceptibility on our time scale is as large as in real 
experiments, the relaxing signal being bigger for x" than 
for x' ■ This explains why the study of x" is usually pre- 
ferred in laboratory experiments. For x" the amount of 



0.25 




Fig. 5. AC susceptibility for L = 100, T = 0.6 and frequencies 
uj=l/P^ 0.02 (circles), 0.01 (squares), 0.002 (diamonds), 0.001 
(crosses). 
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Fig. 6. AC susceptibility for L — 100, T = 0.6 and frequen- 
cies uj = l/P = 0.02, 0.01, 0.002, 0.001 plotted as a function of 
u)t. Following Ref. susceptibilities have been shifted on the 
vertical scale by arbitrary amounts to make them collapse. 
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relaxation is nearly equal to the corresponding stationary 
{ojt — > oo) value. 



5 Conclusions 



In this paper we have made a detailed study of the dynam- 
ical scaling behavior of the correlation functions and AC 
susceptibilities in the aging regime for spin glasses. The 
study has been done applying the scaling behaviors pre- 
ferred by the experimentalists and comparing them with 
the results obtained in numerical simulations of the three 
dimensional Edwards- Anderson model. 

The general conclusion is that there is a full agree- 
ment between the experimental results measured for TRM 
decays and simulations for correlations. This agreement 
must be understood in the following sense. All data col- 
lected for spin glasses is well compatible with a full aging 
scaling scenario with logarithmic corrections. These cor- 
rections are always enhanced in the presence of an exter- 
nal field and they can be misinterpreted as a subaging sce- 
nario. We are not aware of any microscopic model display- 
ing aging with full scaling without logarithmic corrections 
and it is possible (for not to say unavoidable) that such 
corrections are also present in the 3d Edwards- Anderson 
model. This fact explains the small subaging effects mea- 
sured in both experiments and numerical simulations. 

It is also interesting to see how the parameters ob- 
tained in the fits of the decay of the correlation function 
agree with the equivalent parameters for the TRM in ex- 
periments suggesting a universality in the dynamical scal- 
ing of relaxations which is well captured by the Edwards- 
Anderson model. A similar conclusion is obtained also for 
AC experiments which, in general, fulfill a good ojt scal- 
ing law. The collapse of the AC relaxations on a master 
curve, by appropriately shifting them by their stationary 
values shows a nice agreement with experimental results 
and shows how relaxation in disordered systems are de- 
scribed (at least in a very good approximation) by a sin- 
gle timescale (corresponding to the waiting time in TRM 
experiments or to the inverse of the frequency in AC ex- 
periments) . 

Still, the big question must be answered. Why this 
good agreement between simulations and experiments is 
not respected when comparing, at a qualitative level, chao- 
tic and memory effects between the Edwards-Anderson 
model and real spin glasses ? [|3[|2j,^ . It is plausible that 
for time scales short compared to experimental time scales 
some dynamical effects (in particular, chaos and memory) 
are not seen in simulations while the full t/t^ scaling (with 
logarithmic corrections) is already present in this regime. 
Further theoretical, numerical and experimental studies 
are needed to clarify this controversial issue. 
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